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Abstract 



A geometric approach is introduced for understanding the phenomenon of phase synchronization 
in coupled nonlinear systems in the presence of additive noise. We show that the emergence of 
cooperative behaviour through a change of stability via a Hopf bifurcation entails the spontaneous 
appearance of a gauge structure in the system, arising from the evolution of the slow dynamics, but 
induced by the fast variables. The conditions for the oscillators to be synchronised in phase are 
obtained. The role of weak noise appears to be to drive the system towards a more synchronized 
behaviour. Our analysis provides a framework to explain recent experimental observations on noise- 
induced phase synchronization in coupled nonlinear systems. 
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I. INTRODUCTION 

Synchronization phenomena occur abundantly in nature and in day to day life. A few well known examples 
are the observations in coupled systems such as pendulum clocks, radio circuits, swarms of light-emitting 
fireflies, groups of neurons and neuronal ensembles in sensory systems, chemical systems, Josephson 
junctions, cardiorespiratory interactions, etc. Starting from the observation of pendulum clocks by Huygens, 
a vast literature already exists which studies synchronization in coupled nonlinear systems — in systems of 
coupled maps as well as in oscillators and networks [1 & references therein]. In recent times, different kinds of 
synchronization have been classified — mutual synchronization, lag synchronization, phase synchronization 
and complete synchronization [1-3]. 

Many of these studies aim to understand the properties exhibited by the systems once they synchronize 
or exhibit phase-locking. A comprehensive understanding still seems to be lacking when one seeks to 
explain why the systems synchronize. In this paper we introduce a geometric approach in order to address 
this fundamental issue. We aim to understand here the reason for the occurrence of synchronized and 
phase-locked behaviour in coupled nonlinear systems which are subject to weak additive noise. 

We consider a system of n nonlinear oscillators which are coupled. We determine the conditions which the 
dynamically evolving variables of the system must satisfy in order that the various oscillators constituting 
the full system synchronize in phase. Our study at the moment docs not include chaotic systems. We find 
that the presence of weak noise assists in bringing about phase synchronization. 

In section 2, we introduce a geometrical approach to discuss coupled dynamics. We have adapted the 
methods which Wilczek & Shapere developed [4-6] to understand self-propulsion of organisms by shape 
deformations in fluids at low Reynolds number, for discussing the deformations and changes in the orbit 
structure in phase space as the system evolves in time. In section 3 we discuss the dynamics of coupled 
oscillatory systems in the presence of additive noise at the close proximity of a Hopf bifurcation. We show 
the spontaneous emergence of a non-trivial gauge structure for such a system arising from the slow degrees 
of freedom, and induced by the fast variables, and associate it with the geometric approach introduced in 
section 2. In section 4, we obtain the condition required to be satisfied by any two oscillators to exhibit 
phase synchronization & phase locking. Section 5 summarizes the main results of the paper and points out 



3 

directions for future studies. 

Our work is motivated by the need for a theoretical understanding of recent experiments on chemical 
oscillators [7] and numerical simulations [8] which show noise-induced phase synchronization in coupled 
nonlinear systems. (See also [9], [29]). 

II. THE GEOMETRY UNDERLYING THE DYNAMICS OF COUPLED SYSTEMS 

In a series of beautiful papers [4-6], Shapere & Wilczek established a geometric framework to discuss 
the motion of deformable objects in the absence of applied external force. We adopt these methods to 
understand a fundamental phenomenon in nonlinear dynamics, namely phase synchronization in coupled 
dynamical systems. 

In the first section, the main idea underlying the paper is developed, which is based on the approach used 
in [4-6] for a deterministic system. Elaborating along these lines, in the following sections we have studied a 
general system of coupled nonlinear oscillators under the influence of additive Gaussian white noise and we 
find the conditions under which the coupled units within the full system can exhibit phase-locked behaviour 
and phase synchronization. 

We consider a system of n coupled generalised oscillators q(x, t) where the state variables q could in general 
be functions not only of time, but could also depend on a set of additional variables x, say spatial variables 
when there is a metric structure associated with the variables: 

4i = fi(Qi,Q2, ■ ■ ■ ,q n ,IJ>j) , i = 1, ...,n, j = 1,. .. ,p (1) 

Thus qi include also extended systems where the individual elements could mutually influence each other 
through a distance-dependent interaction between the elements. For instance in the case of coupled chemical 
oscillators [8], qi would denote concentrations which have spatial dependence. Another example occurs in 
biological information processing where neurons interact in ensembles. 

We study the simplest case in which, in the absence of couplings, each of the n subsystems admits oscillatory 
solutions for some parameter values fij. Switching on the mutual coupling between these oscillators results 
in the emergence of a collective behaviour. It is then appealing to view the collective behaviour as having 
arisen as a result of some sort of communication between different points in the configuration space. Thus 
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one is led to a geometrical description of the system's dynamical evolution. 

For simplicity, we restrict ourselves in this paper to the case in which the collective dynamics also exhibits 
limit cycle behaviour emerging via one or more Hopf bifurcations. The more general situation which 
includes chaotic orbits for the uncoupled dynamics is not considered here. We define the configuration space 
of the full system as the space of all possible flow lines and closed paths. We consider the situation when 
there is no external driving force, so that the space of all possible contours is the space of oriented contours 
centred at the origin. 

In the absence of mutual couplings, the space of contours consists of oriented closed orbits, each orbit 
inclined at an angle with respect to the other. If we now turn on the mutual couplings between these n 
subsystems gradually, each of the orbits would gradually get deformed, going through a sequence of shape 
changes, and resulting subsequently in a net rotation for it. The problem of interest is to link the dynamical 
variables of the system with the net rotation induced by a change of shape of the orbits for each of the n 
oscillators in phase space. 

The relative orientations of any two contour shapes can be compared by fixing coordinate axes for each. 
Since there exists a degeneracy in the possible choice of axes one can make at each point in the space of 
contour shapes, each set of reference frame we can choose from being isomorphic to E n , a gauge structure 
is therefore induced in this space which facilitates going from one particular choice of axes to another. 

In [4-6], the problem of self-propulsion at low Reynold's number made possible solely through shape 
deformations was discussed. Each choice of reference frame fixed to a shape, which assigned a "standard 
location" in space for each shape, was associated with the motion and location of any arbitrary shape in 
relation to its standard location. We follow their methods closely to discuss deformations of the oriented 
contours in the space of contour shapes. 

Following [5,6], the sequence of oriented contours S(t) can be similarly related to the sequence of the corre- 
sponding chosen reference standard contour shapes So(t) by a rigid displacement 1Z(t) : 



S{t)=Tl(t)S {t) 



(2) 



where in general, an n-dimensional motion 1Z includes both rotations R and a translation I : 

J 




(3) 



where R(t) is an n x n rotation matrix and stands for a sequence of time-dependent rigid motions. The 
contour boundaries are parametrized by the control parameters fa, for each of which the rigid motion 1Z acts 
on the vector [So(/j,), 1] T . The physical contours S(t) are invariant under a local change 

s Q = n[s ]s (4) 

made in the choice of standard contours So- Then the contour shape evolution can be written by combining 
eqn.(4) with eqn.(2) as : 

S(t) = ^(t)n- 1 (S , o(t))5 = K{t)S (t) (5) 



or 



H(t) = n(t)n- 1 (s (t)) (6) 

The temporal change in the sequence of rigid motions can be written as: 

— = Kin' 1 —) = TIA (7 
dt x dt' w 

where A can be identified with the infinitesimal rotation arising from an infinitesimal deformation of So(t). 
Eqn.(7) can be integrated to obtain the full motion for finite t: 



t>/4. \ <nn \<n ft* A{t)dt 
where V stands for the path ordered integral, the Wilson line integral W : 



(8) 



W 21 = Ve£ A(t)dt = l+( A{t)dt + [ [ 

Jti<t<t 2 Jt 1 <t<t'<t 2 J 



A(t)A(t')dtdt' + ... (9) 



in which the matrices are ordered such that the ones occurring at earlier times are on the left. 
It can be seen from eqns.(5),(6) and (7), that A transforms like a gauge potential: 

A = nAn- 1 + n^- (10) 

and the Wilson integral transforms as: 

W21 = fWaif^ 1 (11) 



Shapere and Wilczek exploited the invariance of (9) under rescaling of time, t — > r(t), the measure scaling 
as dt — > rdi, ^4 — > A/f, to rewrite it in a time-independent geometric form. 

This was done [5,6] by defining an abstract vector field A on the tangent space to So. The projection A(t) 
of A at the contour shape So(t) is evaluated in the direction ^ in which the shape is changing: 

A(t)=^ o [£„(*)] (12) 

In terms of these projected vector fields, (8) was rewritten in a time-independent form for a given path and 
independent of the manner in which the path is parametrised in the contour shape space as: 

Rfa) = R(h)vJ*otS A(So)dS ° (13) 

Each of the components Aj[5o] of A coming from each direction in the contour space generates a rigid 
motion and can be defined in terms of a fixed basis of tangent vectors {wi} at So : 

A t [S a } = A m [S a ] (14) 

An infinitesimal deformation s(t) of a contour S (t) can be represented as: 

S {t) = S + s(t) (15) 

where an expansion of s(t) can be made: 

s(t) =^ ai {t)w t (16) 

i 

It was shown in [5,6] that for the particular case Sofa) = Sofa), i,e., for a closed cycle in which the sequence 
of deformations returns the system to the original contour shape in its configuration space, the line integral 
in eqn.(9) becomes the closed Wilson loop which can be simplified to 



W = Ve§ A ^ dt = l + \jY. F ^i 



dt (17) 



where 



The field strength tensor Fij gives the resultant net displacement when a sequence of successive deformations 
is made of So around a closed path and is thus the curvature associated with the gauge potential. 



In the configuration space of contour shapes, the orbit of each of the n subsystems of the full coupled 
systems of oscillators undergoes the shape deformations described above. Because of the mutual couplings, 
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the motion in phase space of any one oscillator coordinate is inseparably linked with that of any other 
phase space point which may be the coordinate of another oscillator. The deformation and motion in the 
configuration space of the various flow lines and closed paths of the entire coupled system can thus be 
viewed as those on the surface of a solid deformable body which is undergoing motion solely due to these 
deformations. 



The full system of n oscillators can be represented by an n-component vector ip l , (i = 1, . . .n) in an 
abstract complex vector space: 

/ \ 

9i 



?2 



(19) 



\q n ) 

A rotation through an angle A with respect to a chosen axis in this internal vector space does not change 
the state of the full system, but just takes one oscillator state ft to another: 



ft -> q t = U (A)ft = e 



it a A° 



(20) 



where t k are k number ofnxri matrices and are representations of the generators of the transformation 
group. Each of the fts represents the state of the ith oscillator at time t. 



There are n independent gauge potentials A Wi corresponding to the n independent internal rotations. Any 
two rotation matrices U(A a ) and U(A b ) do not commute unless A a and A b point in the same direction. 
On application of a common input to the full system, all the different n oscillators respond to it. In this 
case emergence of a collective behaviour is determined by the same gauge potential, although perhaps by 
different amounts or strengths. 

In the following section we would link these gauge potentials with the dynamically relevant variables of the 
coupled system. In section 4 we will attempt to understand the geometrical basis underlying the dynamics 
of phase synchronization between the oscillators. 
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III. COUPLED NONLINEAR OSCILLATORS SUBJECT TO FLUCTUATIONS 

We now consider the system of n coupled nonlinear oscillators % subject to additive Gaussian white noise 
£j in the limit of weak noise: 

?i = /i(?i,?2,---,gVi,M)+& > (« = 1, ••-,"), neR p . (21) 
where the noise correlations are defined as: 

mMt')) = QS lJ s(t-t') (22) 

The eigenvalues of the linear stability matrix of the coupled deterministic system determine the route 
through which the full system moves towards a collective behaviour. A pure imaginary complex conjugate 
pair of eigenvalues at the bifurcation point with the remaining (n — 2) eigenvalues having nonzero real parts 
signals a Hopf bifurcation. The orbit structure near the nonhyperbolic fixed points (qo,^o) OI eqn.(l) is 
determined by the center manifold theorem. When the system described by (22) undergoes a change in 
stability through a Hopf bifurcation, one obtains a p-parameter family of vector fields on a 2-dimensional 
center manifold. 

In this case one observes an emergent common frequency of oscillation for the coupled system. Such a 
situation automatically realises frequency synchronization also since the Hopf oscillator rotates with a 
characterisic frequency. If there are more than one Hopf bifurcations, clearly it indicates more than one 
common frequency of oscillation and one expects to observe a clustering of the various n coupled oscillators 
around these common characteristic frequencies. 

We aim to understand in this paper how phase synchronization results. 

The full system of n oscillators changes stability as the parameter under consideration takes on differ- 
ent values, and at some parameter values undergoes bifurcations. We will study the system in the close 
neighborhood of the bifurcation points where the system exhibits critical behaviour. It is in these regimes 
that the behaviours of the individual oscillators gives way to the collective behaviour of the entire coupled 
system of the n oscillators. We employ center manifold reduction techniques for the system in the presence of 
fluctuations and perform a separation of variables in terms of fast and slow variables as in [10-11], exploiting 
their dynamical evolution on different time scales. A drastic simplification can then be made of the system's 
dynamics and one can write the probability P(qi,t) for the system to be in a certain configuration at time t 



in the weak noise limit as the product: 

P(qut)=p{q f \q a )P(q„t) (23) 

where q s and qf are the slow and the fast variables respectively of the system. The probability P(q s , t) for 
the critical variables occurs on a slow time scale and is non-Gaussian in nature. 

The properties of the fast variables depend upon the nonlinearities in the system. For instance in the case 
when the coupled system exhibits a cusp bifurcation the fast variable could exhibit non-Gaussian fluctuations 
(depending on the specific nonlinear interaction) as it is coupled to the critical variable. It can be shown [10] 
that the joint probability density p(qf\q s ) is confined to a narrow strip peaked about the center manifold. 
We are interested in the case when the coupled system also exhibits self-sustained oscillatory behaviour and 
makes a transition to limit cycle behaviour in the presence of fluctuations. 

It was shown in [10,11], that for a transition via a Hopf bifurcation, p(qf\q s ) has the time-independent 
Gaussian form in the qf variables with width which depends upon the slow variables q s : 

p{qf \q s) = (^M) 1/2 e --te)(w-9/„te)) 2 (24) 

where the center manifold is obtained as a power series in q s : qf = qf Q (q s ). The center manifold theorem 
has been used in [12] for providing with a proof for the stability of the synchronised states. The enslaved 
stable modes are the fast variables which follow the dynamics of the center (critical) modes. We rewrite 
fi(qi,q 2 , ■ ■ .,q n ,tH) as 

, , \ 5F{q 1 ,q 2 ,...,q n ,m) 

fi{qi,q2,---,q n ,Hi) = - x (25) 

oqi 

The Fokker-Planck equation for the full system is: 

Using eqn.(23) we can rewrite this as 

dP(q s ,q f ,t) d dp{q f \q s ) dP{q s ,t) 
dt = J t (P(lf\l S )P(qs,t)) = ^^P(q s ,t)+p(q f \q s )^ r - 

= -(H FPl (qf,q s ) + H F p 2 (q s ,t))P(q s ,qf,t) (27) 

Hence in the close proximity of the bifurcation, the operator Hpp can be written in a separable form, the 
part H F p(q 8 ,t) independent of the fast variables. Here 

-H FP2 (q s ,t))P(q s ,q f ,t) = A(Pfa,t)g) + ^ ( |^ (28) 
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We find it convenient to analyse the coupled dynamics in a path integral framework. We follow the procedure 
of Gozzi [13] to recast the system, eqns.(21),(27) as a path integral, and define: 



* = P{q f ,q a ,t)e F ^)/ 2 



(29) 



so that (26), (28) can be rewritten as 



dt 



(30) 



where 



1 d 2 1 OF 2 ld 2 F 
Wfp - -7:tt-2 + o(-Hr) -7-572" 



(31) 



2 0(tf 8 v % y 4 0<tf 

To enable computation of correlation functions within the path integral formalism, we introduce n external 
sources Ji to probe the full coupled system so that the partition function Z[J] for the system can be written 
as the time ordered path integral 



Z[J] =NT\{ J DqfDq^ie-^S^^^Piqs^piqflq^Siqi-q^e-S 



(32) 



where denote the solution of the system of Langevin equations (21), T denotes time ordering and TV is 
the normalization constant. From cqns.(21) one can write 



5{qi - = 5(qi - /,(gi, ?2, ■■■,q n )- &) 



Sqi 



(33) 



We can rewrite the Jacobian 



<1£l 



5£i 



Sqi 



of the transformation — > qi as 
9fi{qi,q2, • ••,&») 



det 



)s(t-t') 



= cxp{trlna t (%<5(t - i') - d^ 1 -J^-)} 



(34) 



The operator df 1 satisfies the relation 



d t G(t - t') = 5{t - t') 



(35) 



Then we can rewrite eqn.(34) in terms of the Green function in (35) as: 



Sqj 



= exp|tr lnat + In^t-O + Gtf^-O^Ty) } 



(36) 



The system evolves forward in time. Hence 



G(t - t') = 6(t - t') 



(37) 
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Using this and expanding the logarithm in the argument of the exponential, we can simplify eqn.(36) to 



e trln 9t ^>'<K0)aff7y (3g) 



Sqj 

Substituting this back into eqn.(32) and using the mid-point prescription 0(0) = 1/2 of Stratonovich, we 
have 

Z[J] = AfT J] [ DqfDqeD&e-* £ Mt')n(t')dt' J J* a^g-jfc />'(**-/«( 9 i,»,.., 9 „)) a (3g) 

i 

Eqn.(39) can be reduced to 

Z[J] = ATT]l f D g/ D 9s e-*/o J ^*'^(*') d *' e -/o dt '[5^+^ 9 , 2 +^(l^) 2 ] e -^(F W -F(o)) 
i ■* 

= ATTjlf %i)g, e -/o*'[ £ "' + 3 J '( t ')«( t ') e -A( F W- f (0)) (40) 
where we have defined a Fokker-Planck lagrangian 

C FP {q i ,q i ,t) = —q 2 i +—{ — ) +-—— = — q\ + f* - —LL. 41 
4Q 4Q % 2 4Q 2 d^dgi 

which is related to the Fokker-Planck hamiltonian 'Kfp defined in eqns.(30),(31) through a Legendre trans- 
formation: 

l K FP {-K l ,qi,t) = inqi - C FP (qi,q i} t) (42) 

Here iTi are the momenta canonically conjugate to the variables qf. 

SC FP 1 

~5q~ =7Ti= 0* (43) 

so that 

1 i)F 2 1 r) 2 P 1 

= g^ + _(-) + -_ (44) 

We use these relations in eqn.(40) to write the partition function as 

Z[J]=MT\{ f i> 7rj £>g je -/„ t *'[ M "(' r '.«. t )+4 J '(*')*( t ')] (45) 

i J 

From eqns.(27),(29),(30) and (31), we see that in the close proximity of the bifurcation "Kfp can be written 
in a separable form as 

< H F p{-Ki,q i ,t) = , K F p(TTf,q f ;TT s ,q s ,t) = 3i FPl (irf,qf,q s ,t) +Ji FP2 (n s ,q s ,t) (46) 

Thus the corresponding C FP can also be split up as C FP (qi 1 qi, t) = £,f p (qf, qf,t)+ C FP (q s ,q s ,t). Then 
we can write 

Z[J] =MT J ^ / ^ e -/ tdt '[ £ f P («/^/)+^(^^-*)+^^(*') 9 .(*')] (47) 
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An averaging over the fast degrees of freedom enables the partition function to be written in terms of an 
effective Lagrangian as a function of only the slow degrees of freedom. This can be done by first rewriting 
the fast degrees of freedom in action- angle variables (9, 1). The emergence of a non Abelian gauge structure 
can then be seen arising from the evolution of the slow dynamics but induced by the fast variables. After 
tracing the origin of the induced gauge potential to the slow dynamics, we obtain the conditions necessary 
to be satisfied in order for the coupled elements to be synchronised in phase. 

To begin with, we introduce a generating function (qj,I; q s ) which effects the transformation (qj, 717) — > 
(9, 1) to the action angle variables: 

% = 7ri ' aTi = 6i (48) 

S( a \qf,I;q s ) is many-valued and time dependent since the slow variables change with time. 

The phase space structure associated with adiabatic holonomy in classical systems was studied by Gozzi and 
Thacker [14] through Hamiltonian dynamics. We find it useful to employ their methods for our study of cou- 
pled oscillatory systems in a fluctuating environment. Denote Hi(I, q s ,t) = 7ipp 1 (qf(9, 1, q s ),irf(9, 1, q s ), t). 
Using the canonical transformation law, Hi can be expressed in terms of the action-angle variables as: 

dS^Hcif I- a ) 

H x (9, 1, q s (t)) = Hi (I, q„t) + q sl ^ (49) 

Oq sl 

Using the methods of [14] and [15-17], we determine the dynamics of the M critical slow variables q s of the 
system by averaging out the TV fast variables which influence them: 



((Hi)) = — ^ / d N 9H 1 (9,I,q s (t)) 



(27T) 

= TTTiv / d d(H 1 (I,q s )+q sl ^ ) (50) 

(2tt) j oqsi 

where the double angular brackets denote the averaging over all 9: ((f)) = H^y* I d N 9f . 
Since S^ a \qf,I;q s ) is multi-valued, the single-valued function 

C(M,«a) =S^(q f (9,I,q s ),7r s (9,I iqs ),q s ) , (0 < 9 < 2ir) . (51) 

is introduced [18]. We have 

d^ = os^ + 7r dQfi (52) 

dq si dq sl U dq sl 
Hence this can be substituted into eqn.(50) to obtain 

{{H,)) = Hi(I, q s ) + q sl ((P~ TT^^i)) (53) 

Oq sl Oq sl 
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The total Hamiltonian of the system is given, after performing the angle averages by: 

H av (I,n s ,q s ) = ((5Ci(g/,7r/;g a ) + { M 2 (q s ,it s ))) = {{Hi(6,I;q s ) + W 2 (q s ,ir s ))) 

= S(I M ) + 4sd(^-^)) (54) 

where we have let 

H(I,n s ,q s ) = H 1 (I,q s )+X 2 (q s ,ir s ) (55) 
The Gibbs partition function in eqn.(45) can be rewritten in terms of the fast and slow variables as 
Z[J] = NT j D-K f D-K s Dq f Dq s e~ So ^'l^fPiiij ^f-gs)+KFP 2 (<]s^ s )+^(J4t')q s (t')+j J (t')q J (t'))] 

= NT j D 7 r s ^ s Z?/^ e -/o d *'^( / ^'^)-^^^^ + ^( J ^ t ')^( t ') +J ^*')^ t '))] (56) 

Performing the integration and simplifying the resulting expression, we get 

Z[J]^NT J Dn s Dq a DIe-S° dt ' {Hav{I ^ s ' qs)+J z Jsqs) (57) 

We use the Magnus expansion [19,20,21] for expanding the time-ordered integral, which gives the final state 
properties in terms of integrals over the initial state ones, to rewrite eqn.(57) as: 

Z[J] w N J DTT s Dq s DI «xp{-^ dt' [h av + s q s + ^[H av {t'), J H^t^dh] 

+ ^[H av (t'),J [H av (t 2 ),J H av {t 1 )dt 1 ]dt2] + ^[[H av (t'),J H av (t 2 )dt 2 ], J H„„(ti)dti] 

+ •••)} (58) 

The variation of parameter(s) in in time of the system brings about change in its stability. The commutator 
terms in the Magnus expansion hence arise on account of this parametric time dependence of the Hamilto- 
nian: \H av {t),H av (t{)] = [H av (m(t')), H^iU))}. 

The necessity of the time ordering is also motivated by the work of ref. [22] who have shown that the order of 
arrival of signals at an oscillator in a network of pulse-coupled oscillators is crucial in determining changes 
in its phase. 

Retaining terms only upto the first commutator in the expansion and substituting for H av from eqn.(54), 
we obtain after some simplifications: 

Z\J\ « M ( D7r s Dq s DI cxp { - / dt'(H(I,w s ,q s ) + g s ,((^-» - q Sl {(« j + ±-J s q s 

J L Jo a( ls i oq sl y 
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- hH{I,ir s ,q s ,t'), I H(I,Tr s ,q.,t 1 )dt 1 ] + hH(I,Tr s ,q a ,l/), f dt iqSm ((n fk ^))] 

Z Jo L Jo °1s m 

+ \[is m {^f k ^)), f H{I,Tt„ qa M)dti\-\[q.A{*U^)), f dt iqsm ((n fk ^))}\ (59) 

z °1sm Jo z °1sl Jo °1s m J 

From here we can define an effective Hamiltonian H e ffi 

F e// = ^(7,^, gs ) + ^((^))-^((^.^))-i fe (( 7r/ .^i)), f dtiqSm{{7Tfk ^I±))] (60) 

oq sl Oq sl Z oq sl J a a Qs m 

From a Hamiltonian variational principle, it was shown in [14] from simple arguments that the averaged fast 
motion induces an effective gauge field which acts on the slow variables. We follow these arguments closely 
for the coupled system subject to fluctuations near the instability. The variational principle gives: 



SS, 



:ff = 5 j dt[ir sl q sl - H ef f(I,ir s ,q s )] = 
Jo 

= 5 ^ T di{[^ + (< 7 r /i |g))]^-5(/, 7 r s , gs ) + i[g S! (( 7 r /i |^)),^ d tl q Sm ((n fk ^-))}} (61) 

The term having the single-valued function ( vanishes since it is a total time derivative. 
Varying S e ff with respect to n s and q s , keeping the end-points fixed gives: 

-a^-'"]} (62) 

We define as in [14], the quantity in angular brackets as 



^JuT ))=Ai (63) 



Then SS e ff = leads to 

dfi 



Qs, = 



dir s 

dH (dA m dAi I 

As in [14] we can identify Ai with a gauge potential, and a curvature tensor Ti m can be defined as 

dA m dAi 1 
fim = -^ o V-[A u A m ] 65 



so that the momenta in (64) can be rewritten as 



"o |--^mQ 66 



15 

The commutator terms in the momenta and curvature tensor were absent in rcf.[14] since the Magnus 
expansion for the time ordered integral was not used there. 



It should be noted that these commutator terms in eqns.(61), (62),(64)-(66) arising from terms such 
as [H (n) , H (/j,')] in eqn.(58) could in general impart a non-flat nature to the connection and generate a 
curvature which is non-trivial by effectively generating new parameters which were not present in the original 
Hamiltonian. (Sec for example [23,24] where in the context of examining the connection between classical 
and quantum anholonomy for some interesting systems (in particular the displaced harmonic oscillator 
[25]), it was shown that for time- varying Hamiltonians, the original Hamiltonian must be embedded into a 
larger class for locating the effective parameter space where the Berry phase two- form has singularities). 
The appearance of a nontrivial gauge structure in general dynamical systems, including classical systems, 
due to a slow variation of the parameters was also explicitly demonstrated in the seminal work in [26] . 



Eqn.(64) shows that the curvature tensor T\ m exerts a velocity dependent force on the slow variables. In order 
to write canonical equations of motion, one has to therefore introduce modified Poisson bracket relations in 
the slow-variable space: 

, , df dg df dg , r df dg 

if {q s ,n s ), g{q s , Tr s )\ = {- — — - — } - T im - — 67 

OK si Oq sl dq sl dTTsl dlt sl dir sm 

Thus the gauge potential coupled to the slow variables is induced by the fast degrees of freedom as is 
evident from (63), the spontaneous appearance of the gauge symmetry being associated with the phase 
degrees of freedom of the center modes. 

The emergence of a gauge structure for the system follows from the crucial property of separability of the 
variables as slow and fast ones evolving at different time scales, which results from the slaving principle 
for the stable modes near the bifurcation in a noisy system. This leads to the motion and deformation 
of the closed orbits in the configuration space. The rotational symmetries of the sequence of successive 
deformations of each orbit brought about the gauge potential discussed in Section 2. The analysis above in 
the current section shows that this can be related to the dynamically evolving variables of the coupled system. 



Having made the correspondence of the gauge potential Ai and the curvature tensor T\ m with the 
dynamics of the actual coupled system through the fast and slow variables, we proceed to examine under 
what conditions phase locked behaviour and full synchronization would occur in a coupled system. 
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IV. CONDITION FOR SYNCHRONIZATION BETWEEN THE COUPLED OSCILLATORS 

At any instant of time, the phase difference between two oscillators q\ and q 2 located at two different 
points x and y in the configuration space can be found from their inner product: 

_ (92(2/), 91(2/)) 
v ~ \(q2(y)Mv))\ 

= [ d d yd d X Tr(p(ei: A ^ tads n) j g f !' gl( /l (68) 

J v ' I (92(2/), 91(2/)) I 

where 9 y denotes the angle between the oscillators qi(x) and q2(y) in configuration space, measured at the 
coordinate y. The path-ordered Wilson line integral appears in the equation above since qi(x) must be 
parallelly transported to the coordinate point y in order to compare it with q 2 located at y. P denotes the 
path-ordering. Since the qiS are related to each other through a gauge transformation in the n-dimensional 
configuration space, this can be rewritten using eqn.(20) as 

(qi(y),qi(x)) 



cos 



e y = ( rfVWe-^ p{j: ' , gl y >gl * , 

J v ' I (92(2/), 91(2/)) 

= ( ^^Tr(e-^ A V(e/:^W ta ^)P(ei';^ (rt ^^) N ) 

V ' 1(92(2/), qi (y))\ 

rf d ^Tr(e-^ A3 P( e r A " (s)ta ^)P(er^ (p)t7 ^)) ^ -i (69) 

V 7 I (92(2/), 9i (2/)) I 



since (qi(y), qi (2/)) = 1. If the angle between gi and 52 remains constant for all times, then the oscillators 
<7i and 92 would be phase-locked; if the angle between them is vanishing for all time, the oscillators would 
be fully synchronized in phase with each other. 

We would like to determine the conditions under which the phases of any two oscillators in a coupled 
nonlinear system would be locked and fully synchronised. Since each ^ is an oscillator, each undergoes 
periodic dynamics in the configuration space. 

Let qi(y) after being parallelly transported from coordinate x, now return to the point x during the course 
of its temporal evolution. We denote the state of this oscillator after completing one orbit and returning to 
x by q[(x). By the time q\ completes this orbit, q 2 would have evolved to another point z. Hence we would 
now like to calculate the angle between q% (z) and q[{x). 
We have: 



9i 



[(x) = P{e§ A ^ s)tads ») qi {x) (70) 
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q2(z) can be parallelly transported to x to compare it with q[(x): 

q2(x)=P{J: A ^ d ^) q2 (z) (71) 

Then the angle between qi and q 2 at x can be calculated: 
0fc(z),3i(*)) 



COS tfj; 



IM*),^*))! 

^ V y I {Q2{x),q' 1 (x)) | 

/^z^Trf^le/r^^^^^e^^^^Ce/^W^^))-^^ 1 ^ 
J V y |(«2(ar),g((x))| 

| d d zd d xTr ( (P( e r A * ds » y**^ P{Jl A ° (,) "'»)) T P(e^ A » )) • 



(91 (z), 9i (a;)) 



(72) 



iM^lC*))! 

The change in the angle between q\ and q 2 during a time interval t can be found using eqns.(69) and (72) 
and by simplifying the resulting expression, to lowest order in A to be 



cos 8 y — cos 6 X — — 



/ ( f x f Paa fx 
d d xd d y[ - ik a F^ + F% J A^(s)ds fl + (S aa H —A^)F^ J A^(s)ds^ + ...} 

(73) 



In arriving at this expression we have used the relation for the Wilson loop integral 

P(J = e F,~ (74) 

in which the t a (introduced earlier in eqn.(20)) are generators of the Lie algebra 

[t a ,t b ] = ie abc t c (75) 

and 

t a F; v = d^Al - - i[t a A«,t b Al] (76) 

is the curvature tensor of the complex abstract vector space. Also we have used the matrix identity 

e A e B = e A+B+±[A,B]+... (7?) 

and the trace relations: 

tr(i£) = 

tr(t»tj) = C(p)S ab (78) 



where C(p) is a constant for the representation p. As the system we are considering is subject to fluctuations 
and is not deterministic, the quantity which is actually of interest to us is the noise average {cos9 y — cos9 x ) 
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of the phase difference between qi and qi: 



(cos9 y -cos6 x ) = p( / d d xd d y{ 




+ 




f A*(s)ds, + (6™ + ^-hP)Ff. f A«(s)d s , + ...}) 




(79) 



Using the Gauss-Bonnet theorem, we see that the first integral on the right hand side of this equation gives a 
topological invariant, the Eulcr characteristic xe of the surface 5* over which the integration is performed: 
Is Fij — Xe- For the situation in which there is perfect phase synchronization between any two oscillators 
in the system, this constant term on the right-hand side of eqn.(79) should vanish and the other terms in the 
equation must also vanish. The two-torus T 2 is a well known example of a topological space with vanishing 
Euler characteristic. The limit cycles of the coupled system are therefore constrained to remain on T 2 as 
they synchronize in phase. 

For the oscillators qi and q 2 to exhibit phase-locked behaviour, we observe that we must have, to lowest 
order in A, 



While the explicit value of the left hand side of (80) would vary from one set of coupled systems to another, 
it is interesting to note that in all cases external noise seems to play a role in bringing about the phase 
synchronization. This can be seen as follows. The noise averages of and are calculated by solving 
the coupled Langevin equations in (21) using eqns.(63) and (65). From eqns.(22),(48) and (63) we see that 
the lowest order terms of (A^) and {F^ v ) such as ( ((tt/J) } would not contribute so that at least to 
this order, the left hand side of eqn.(80) is brought very close to zero taking the system towards synchrony, 
whereas in the case when no fluctuations are present these terms would be non-zero. 

The role of noise in bringing about phase synchrony can also be understood in the following way. Consider 
the analysis by Ermentrout [27] of two weakly coupled oscillators: 



Zi E R Ni , Fi are continuous and differentiable, Gi are continuous, and each uncoupled system — Fi(Zi) 
admits a unique, globally stable periodic solution. It was shown in [27] that the coupled state admits a 
parameter regime in which n : m phase-locking occurs between the two oscillators, after n cycles of oscillator 
1 and m cycles of oscillator 2. Using a multiple-scale perturbation technique, introducing slow r and fast 
s time variables: ^ 

= ^Ts + K 37 > this + ^-dimensional system was reduced to a one-dimensional 




(80) 



— ^Tf- = Fi(Zi) + nGi(Zi,Zj) , i,j = l, 



2, % + 3- 



(81) 
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evolution equation on a slow time-scale for the phase difference $ between the two oscillators. This was 
shown to have the form: 

— = H(<P) (82) 

where H (<£) = + 2n), and the phase shifts <£>s vary slowly in the direction of the flow of the limit cycle 
which is formed due to the coupling of the oscillators. Further it was shown that the phase-locked solution 
to the coupled system corresponds to the fixed points of eqn.(82), H being identified with the Poincare map 
for the flow of the full system. 

For a nonzero value of the noise strength, the sharp transition to the critical point is replaced by a 
bifurcation region, and hence the time spent by the unstable modes near the bifurcation (critical slowing 
down of the deterministic system at the bifurcation) is much longer in the stochastic case. 
It is known [28] that the slower the system moves along any part of the limit cycle, the larger is its statistical 
weight in that part of the limit cycle. Hence, from eqn.(82) and the result of [27] mentioned above, the 
phase-locked solutions of the coupled oscillators are statistically favoured. Moreover, since in the stochastic 
system these (unstable) slowly varying phase differences show an increase in the relaxation time as the 
instability is approached, as compared to the deterministic case, phase-locking and synchronous solutions 
have a larger statistical weight in the presence of the weak noise. 

We have considered the case in which the coupled system exhibits limit cycle behaviour. The formation of 
a limit cycle involves symmetry-breaking, permitting the existence of both stationary and time-dependent 
probability densities (for finite and infinite system volumes respectively) . The probability peaks for the time 
dependent densities rotate along the limit cycle, while the time-independent densities are crater shaped. A 
particular phase is associated with every given realization of a limit cycle and a choice made corresponds to 
breakdown of gauge symmetry. The emergent gauge structure associated with the phase degrees of freedom 
of the center modes in the vicinity of the bifurcation, enables us to introduce the geometrical quantities A\ 
and Fi m , and to obtain the condition for phase-locking (eqn.(80)) in terms of these quantities. 

Our analysis was made possible only because of the presence of weak noise which, as we showed in Section 
3, plays the crucial part of enabling separability of the variables near the bifurcation into slow and fast ones 
evolving at different time scales. 

Measurable phase differences between oscillators in the presence of the (weak) Gaussian white noise are 
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noise-averaged quantities, and from eqn.(80), are determined largely by the fluctuations in the gauge 
potential. Since these average out to zero, we conclude that the presence of (weak) Gaussian white noise 
always enhances phase synchrony. 

Recent experimental observations by Fujii et al [7] of two chemical oscillators separated by some distance 
in the light-sensitive Belousov-Zhabotinsky reaction show self-synchronization of phase and frequency by 
application of noise. They observed spontaneous synchronization for small separation distances in the ab- 
sence of noise and demonstrated the existence of an optimum noise intensity for the self synchronization 
phenomenon. Phase synchronization in coupled non-identical FitzHugh-Nagumo neurons subject to inde- 
pendent external noise was also demonstrated through numerical simulations in [8] . Noise-induced phase and 
frequency synchronization was also demonstrated recently in stochastic oscillatory systems both analytically 
and with numerical simulations [9]. (see also [29]). 

Our analysis provides a framework to understand these findings and opens avenues for deeper studies 
relating cooperative phenomena in coupled nonlinear stochastic systems, with the underlying rich geomet- 
rical structure of the phase space generated by the complex dynamics, and with associated mathematical 
invariants which govern the system's asymptotic behaviour. 



V. CONCLUSION 



We have introduced a geometrical approach aiming to understand phase synchronization among coupled 
nonlinear oscillators subject to additive noise. We have considered the specific scenario when the collective 
dynamics of all the oscillators also exhibits limit cycle behaviour arising via one or more Hopf bifurcations, 
consequently implying the occurrence of frequency synchronization. We demonstrate the emergence of a non 
Abclian gauge structure arising from the evolution of the slow dynamics but induced by the fast degrees of 
freedom. The condition required to be satisfied in order for phase locking and phase synchronization to be 
exhibited is obtained in terms of characteristic invariants of the surface generated by the dynamics of the sys- 
tem. We find that weak noise helps in bringing about phase synchronization. This provides an explanation of 
recent experimental observations and numerical simulations of noise-induced phase synchronization [7,8] (see 
also [9,29]). Our work also motivates further studies of the internal structure and geometry of synchroniza- 
tion defects in spiral waves in oscillatory media which have been areas of keen interest in recent times [30,31]. 
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